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1 Introduction 

Hofava has recently proposed a theory of gravity which is power-counting renormahzable and 
possibly unitary [1]. This is achieved by assigning different scaling dimensions to time and 
space, at the cost of breaking Lorentz and diffeomorphism invariance at high energies. In four- 
dimensional Hofava gravity, one usually considers an anisotropic scaling exponent oi z = 3, such 
that t ~ x^. This means that the Lagrangian should include higher spatial derivative terms up 
to sixth order. One may also consider couplings to scalar [21 [3l HI [5] and vector [21 |5l El [7] fields 
with higher derivatives terms. There is an ongoing discussion regarding possible flaws of the 
theory. We refer to the latest reviews on these subjects [5], and references therein. 

The study of Black hole solutions and thermodynamics in general relativity with higher 
curvature corrections has been pioneered by [9]. There, large black hole solutions were obtained 
by treating the higher curvature terms as a small perturbation. Black hole solutions in Hofava 
gravity and their thermodynamics have been discussed in [lOl [HI [121 [131 [IH [IHl [161 [H]- ^^ 
particular, [HI [TB] discusses charged black holes in Hofava gravity coupled to the usual Maxwell 
term. In this work we add higher derivative magnetic terms to the mix, and discuss black hole 
solutions and thermodynamics. 

The paper is organized as follows: In section 2, we discuss Hofava-Maxwell theory with 
higher derivative magnetic terms. In section 3, we discuss static spherically symmetric black 
hole solutions in the low-energy approximation. In section 3.1, we calculate the horizon locations 
and temperatures for near-extremal asymptotically fiat black holes. In section 3.2, we do the 
same for near-cold asymptotically (anti-)de Sitter black holes. In section 4, we discuss a detailed 
balanced version of the theory, and find non-perturbative solutions. 



2 Hofava-Maxwell Theory 

We use the ADM formalism where the metric is parameterized as 

ds'^ = -N^df + Qijidx' + N'dt){dx^ + N^dt) , (1) 

where i,j run over spatial coordinates. The extrinsic curvature is given by 

K^J = ^{g^j - V.iV, - V,N,) , K^K\, (2) 

where dot denotes a derivative with respect to t. The four- dimensional Ricci scalar may be 
decomposed (see e.g. [18]) as 

'■^^R = R + KijK'^ -K^ + (covariant derivative) , (3) 



where R is the three-dimensional Ricci scalar. The last term contributes a total derivative in the 
Einstein-Hilbert Lagrangian, and may be dropped in this case. Note however, that this will not 
be true in general for higher curvature terms. Thus, Hofava gravity, which is of higher curvature 
in R, will differ from the more common ^^^R higher curvature theories, even for static solutions. 
The Cotton tensor is given by 

ikl / -| 



C., = ^V. \^R\ - -R6jj , (4) 

where e^^^ = 1, g = det{gij), and Rij is the three-dimensional Ricci tensor. 

Hofava's "detailed balance" Lagrangian density [Ij is given by: Cuofava = ^o + ^i, with 



>Co = -,{K,,K^'-XK') + ^rY^MwR-^^'w) 



N^ " k2' ^ ' 8(1 -3A) 

' -C, = ^-^il^i?^ - ^^R^^R^ + ^'—R,V,R^, - -^r,C- , (5) 



N^ 32(1 -3A) 8 ' 2w^^ ' " 2w^ 

where k, A, /i. Aw and w are constant parameters. Comparing Cq to the Einstein-Hilbert 
Lagrangian, these parameters are related to the speed of light, Newton's constant and the cos- 
mological constant by 



We want this theory to reduce to general relativity in the low-energy limit. We therefore first 
set A = 1. Choosing units so that c = 1, shows the coefficients of the R^, RijR^^ terms in £i 
to be inversely proportional to Aw, with no additional free parameters. In order to allow both 
the cosmological constant term and the higher curvature terms to be small, we must deform 
the theory to break the relation between these terms |1], |10]. As long as we do not change the 
form of the six- derivative CijC^^ term, this is considered a soft violation of the detailed balance 
condition. We will write this as: Cdejormed = ^o + ^clqAwCi, where Oq is a constant parameter, 
and e > will be used later as a small expansion parameter. In this parameterization Av^^ factors 
out, leaving the now independent usual cosmological constant term (see e.g. (fT5|) ). 

We now consider electromagnetic interactions. Maxwell's Lagrangian density in the ADM 
decomposition takes the form: 

^£2 = ^,9''{Fu - N'^Fu.){F,, - N'F,,) - F,,F^^ , (7) 

where 7 = IGnGc^^ fixes the normalization in relation to the gravity part, and with the field 
strength given by 

Fu = Ai - diAt , Fi, = d,Aj - d,Ai , (8) 



where [At, Ai) is the vector potentiaL We will add minimally coupled higher derivative magnetic 
terms as considered in O [6] : 
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^, -£3 = -aiVkF^^VT^' - a2VkViF,,V''VF'' , (9) 

where ai, 02 are constant parameters. [6] discusses constraints on the magnitude of 02. Note that 
contrary to the case of O |6] which consider cosmological backgrounds, terms related by integra- 
tion by parts to the RHS of iQ are relevant since A^ is space-dependant (in the non-projectable 
case). Nevertheless such terms are not considered here. Note also that in the spherically sym- 
metric case considered later, VjF*-' = and AF*-' = ^e^^^ e^"^"^ R^iF^n- Therefore, we have not 
considered here terms containing such derivatives. In section 4, we consider a detailed balanced 
version of Hofava-Maxwell theory, which does include such terms. 

We take our overall Lagrangian to be 

£ = £0 + eaoAH/£i + £2 + ££3 • (10) 

Finally, we make the large black hole (or low-energy) approximation, in the spirit of [9]. We 
require that the gravitational and magnetic higher derivative terms are small compared to R and 
FijF^^ , respectively. In our analysis this will usually mean 

e|«o|<r^or> e|ai|/<r^„^, e |a2|/ < r)j„, , (11) 

where Vhor is the location of the horizon of interest, and p is the magnetic charge. However, for 
asymptotically anti-de Sitter solutions with a large (negative) cosmological constant: \Aw\ rf^^ ^ 
1, the conditions instead read 

e\aoAw\<^l , e|ai|/<r^„^, e |a2Ai^|/ < r^^r • (12) 

This may be seen from fl20|) . We will construct perturbative solutions of this theory to first order 
in e. Our calculations were done using Maple with GRTensor. 



3 Static Spherically Symmetric Black Holes 

We consider static spherically symmetric solutions, where the metric takes the form: 

ds^ = ( - N(r? + NJrff(r)]dt^ + 2NJr)dtdr + -^dr^ + r\de^ + sin^9d(f)^) , (13) 

V J f{r) 

where we have retained the shift function A^,. as discussed in [T3| [H] IT5] . The vector potential 
takes the form: 

At = A{r), Ar = Ag = 0, A^ = -pcos9 . (14) 






Substituting the above ansatz into the Lagrangian density we get 

a/7 

yJ-r^Creduced oc SA^yr^ + 2/V^ + 2fr' - 2r' + eao{2f'fr' - 2f'r' - fr^ + 2fr^ - r') + 

- ^A'^r^ + 2pV^ + Ueaip^fr^ + 3ea2P^(/'V^ - Uf'fr + 60/^) + 

+ ^(^KNrfr' + N'j'r' + N'jr') , (15) 

where prime denotes a derivative with respect to r, and we have omitted an overall factor inde- 
pendent or r. Varying Creduced with respect to A^, Nr, f and A, yields the equations of motion, 
the latter giving 

where q is the electric charge. 

The variation by A^ gives 

A, 

N \N f 

Let us first consider the case Nj. ^ 0. The motivation here is to obtain solutions to the projectable 
version of Hof ava gravity, where A^ = N{t) [H [131 [IB] • In our static case, A^ must then be constant, 
implying that also / is constant due to (ITTl) . For ai = a2 = 0, uncharged and charged solutions 
of Nr were discussed in [13] and [16], respectively^ For general ai ^ or a2 ^ 0, we find 
there are no magnetically charged solutions to the equations of motion. ^E\ discusses black hole 
thermodynamics when A^^ 7^ 0. From here on we assume A^^ = 0. 

The zeroth order solution is the Reissner-Nordstrom-(anti-)de Sitter black hole: 

Ar. = ;„._^,.+ l_H^ + ^, (18) 

where mo is the mass at zeroth order (see however ^7\). For ai = a2 = 0, uncharged and charged 
non-perturbative solutions were obtained in [TO] and [UllIE], respectivelyo For the general case, 
we will look for solutions to first order in e, in the form: 

/ = /o(l + eyp) , N' = /(I + er]) . (19) 

Note that in the expression for A^^ we use / rather than /o, and truncate above order e. The 



^See e.g. (j55|) with (3 — 0, assuming ea^Aw — 1- 
^See e.g. (|56l) with /? = 0, assuming eaoAw = 1. 



solution reads 



V 



105ao (A^r^2 - AKw{q^ + P^V + IQmlr^ - 16mo(g2 + p^)r^ + A{q^ + p^) V^) + 

+ 168aip^ (ihKwr^ - lOr^ + 15mor^ - 6(g2 + p'^)r^) + 

+ 8a2p^ ( - 1575A^r^ + 2520Aw'r^ - 4095Avi^mor^ + UGAAwiq^ + p^^ - 1890r^ + 

+ 6930mor^ - GSTOm^r^ - 3240(g^ + p^)r^ + 6300mo(g^ + /)r + 

-1540(g2 + p2)2) + 



and 



+ leSOmn^r^ 



r] = 3p 



420(Awr^^ - 2r^° + 4mor^ - 2(g2 + p^)r^] 



—ai + 9a2Ai^ 



10a2 9Qa2mQ 



4a2(g^ + P^ 



(20) 



(21) 



where 5 is an integration constant, and we have assumed the e-corrections vanish at infinity. Since 
in general rj y^ 0, N"^ is not proportional to /. Thus there is no transformation of coordinates that 
would yield a corresponding projectable solution [13], as indeed expected in light of the above 
discussion. 

Looking at A^(r)^, one sees that the 5-term is just an explicit e-correction to the mass: m = 
777-0 (1 + ^S). We assume the horizon locations of the black hole are given by the roots of gu = 
N{rY = 0, as in the relativistic case (see ^14j for a discussion of the notion of horizon in Hofava 
gravity). One may set 6 to make f{r) and r]{r) nonsingular on the unperturbated black hole 
event horizon, as done in [9]. This would mean the event horizon location does not get an e- 
correction. However, the mass then has a fixed e-correction, and this depends nontrivially on the 
zeroth order mass rriQ. Also, other horizons of the solution do get e-corrections. We will take 
the alternate route of leaving 6 free, and finding all horizon locations perturbatively, designating: 
''^hor = ^o(l + £p)) with e \p\ <^ 1. Note that in the near-horizon region, /(r) and N'^{r) now get 
an additional e-correction due to p in /o(^), and it is the overall correction that should be small. 

The Hawking temperature of our static spherically symmetric black hole (with A^,, = 0, 
otherwise see [U]) is given by 



T 



4vrV-5'tt^r 



27r 



^ '^hor 



(22) 



"^ '^hor 



Note that for asymptotically de Sitter black holes, the temperature may need to include the 
Bousso-Hawking normalization factor [19j, not considered here. We do not calculate the Bekenstein- 
Hawking area law entropy, as it is expected that the entropy would be corrected by the higher 



derivative terms. In fully diffeomorphism invariant theories, this is computed by the Wald en- 
tropy [20]. However, as mentioned also in |14j, it is not clear how to do this in our diffeomorphism 
breaking theory. In the following we will analyze the horizon locations and temperatures for some 
simplified cases of black holes. 

3.1 Asymptotically Flat Black Holes 

We first consider asymptotically flat black holes with Ky/ = (holding c finite). The zeroth order 
inner and outer horizons are given by 



Tq" = mo- yml - q^ - p^ , r^ = mo + \J ml - q^ - p^ . (23) 

The expressions for the first order corrections are lengthy, and to be concise we will consider only 
the extremal and near-extremal cases. The extremal limit is determined by: r™ = r°"*, which is 
satisfied at zeroth order by 

r^^*'^ = ml'''' = ^/^fTf . (24) 

At this extremum point, the first order correction to the horizon location will not affect N{r)^ 
to first order in e. In order to satisfy A^(r)^ = 0, one must set 

.extr ^ -105ao(g^ + p^Y + 42ai(g^ + p^)p^ + 20a2P^ 

420(g2+p2)3 ' ^ ^) 

correcting the extremal mass. It is conjectured that the mass-charge ratio of extremal black holes 
is decreased by higher curvature corrections [21j, implying: 5^^*'' < 0. The first order correction 
to the extremal horizon may then be got either by taking the limit of the non-extremal case, or 
by satisfying: (A^(r)^) = 0, that is, having the temperature vanish. It reads 

Note that the extremal horizon does not depend on oi, a2- 

In the near-extremal approximation we parameterize the mass as 



m = m^^*Vl + A2 , m"^*'' = m[;^*"(l + e<5"^*") , (27) 

and expand all expressions to leading order in A ^ 1. The horizon locations satisfying: A^(r)^ = 
+ O(A^), are given at zeroth order in e by 

r^ = V^T^il - A) + 0(A2) , r°"* = ^/^T^il + A) + ©(A^) , (28) 



and at first order in e by 

^n ^ -105ao(g^+p^)^-336ai(g^+pVA-100a2P^A ^ 

^ 420(g2+p2)3 +^IAJ, 

out _ -105ao(g^ + p^f + 336ai(g^ + p^^ A + lOOa^P^A , ^, ^ ,, 
^ " 420(^2+^2)3 +cy(Aj. [ZJ) 

The Hawking temperature is 

"^ (^ I .525ao(,2+p2)2-294a,(,2+pV-80a.p2^ ^ ^^,, ^ ,^ ^3^^ 



27rv/^^+^ V 420(^2 + ^2 

3.2 Asymptotically (Anti-)de Sitter Black Holes 

We now consider asymptotically (anti-)de Sitter black holes with Aw 7^ 0. In the de Sitter case, 
having Aw > (holding c positive by analytic continuation), there exists a cosmological horizon 
denoted: r™'^"^. Let us recall the classification of the limiting cases of Reissner-Nordstrom-de 
Sitter black holes [22l[T9l: 



(a) lukewarm: r'"" < r°^^ < r""'"" and T'"'* = T'^"^™. 

(b) charged Nariai: r™ < r""* = r'^"*™. 

(c) cold' r*" = r°"* < r'^"^^ 

(d) ultracold: r*" = r°"* = r'^°'^™. 

In the anti-de Sitter case. Aw < 0, and the only relevant limit is the cold black hole: r*" = r°"*. 
Together with the fiat case, these limiting cases (and their "near" versions) have a concise zeroth 
order form for the horizon location. We will first consider the cold and ultracold limits. We will 
then extend to the near-cold black hole: r*" < r°"*. For the latter to be concise also at first order 
in e, we will consider the near-fiat approximation: \Aw\ {q^ + P^) ^ 1- 

The zeroth order cold black hole metric function is of fourth degree in r, and has a double 
zero. It may be written as 

W)' = -^(^ - ^)(- - Bf{r - C) , (31) 

where A, 5, C are constants. Equating fl3T|) with flTSj) and comparing powers of r, one gets 

A + 2B + C = , 
SAwB"^ + 2AwBC + AwC^ -2 = 0, 
AwB^ - B + mo = , 
3AwB^ - 2B^ + 2(g2 + p^) = . (32) 



This gives 



and 



^cold ^ ^ 



2(g2+p2 



I - 6Aw{q^ + P^) + I ' 



m. 



cold 



|(2 + v/l-6AH^(g2+p2; 



where for Aw > we have assumed B < C, and in that case we have also 



C = -B 



U + 2^1-6Aw{q^+P^) 



3A 



w 



(33) 
(34) 

(35) 



As in the extremal flat case, to satisfy A^(r)^ = to first order in e, one must correct the cold 
black hole mass by setting 



6 



cold 



70ao(g' +/)'(- 7 + 12AH/(g2 +/) _ 2y/l - 6Aw{q^ + p^)^ + 

+ 126ai(g2 + p'')p^ (^1 - 18Aw{q^ + p') + y^l - QAw{q^ + P^)) + 
+ 3a2p' (20 - 318AH'(g' + p') + 2844A^(g2 + p^ - 23y/l-QAw{q 

+ 43{l-6Aw{q'+p')f' 



2+p2) + 



l + y/l-QAw{q^+P^ 



1680(g2 + p2)3 2 + ^1 - QAwiq^ + P^ 



The first order correction to the cold horizon location reads 



^coid „ l + ^l-6Awiq'+p^) 

8(g2 + p2)^(l_6AH.(g2 + p2 



P 



-oo 



(36) 



(37) 



Here again the cold horizon does not depend on oi, 02- The full expressions for the cosmological 
horizon location and temperature are more cumbersome, and will be given later in the near-flat 
limit. 



The ultracold limit is determined by: r^"'*^ = r'^"*™^ which is satisfied at zeroth order by 

1 



A, 



uc 



6(g^ + p"^ 



{31 



and 



MC 



-m. 



uc 




2(g2+p2 



(39) 



where we have denoted: Aw = Ao(l + ev), and we expand in e before taking the ultracold limit. 
Note that the first order corrections to the mass and horizon location, now also get a contribution 



from their respective zeroth order terms due to v. Thus the e-correction in the expansion of fl33|) . 
serves to cancel out the divergency of fl37|) in the hmit f l38|) . The overall first order terms are 
given by 



^uc _ '^0 



2(g2 + p2) ' 



jjc _ -280ao(g' + p'f - 126ai(g' + p'^ + GQaap' 



,2-\3 



1680(g2+p 
In the near-cold approximation we parameterize the mass as 



m = uf"^'^ 



V1 + A2 , nf°^^ = mf'^{l + e6^°^'^) , (41) 



i2 



and expand all expressions to leading order in A <^ 1. The horizon locations satisfying: N{r 
+ 0(A'^), are given at zeroth order in e by 

r^" = 6(l-a,A)+0(A2), r"''' = b{l+abA)+0{A^) , r^"'"" = cil+acA^)+0{A^) , (42) 

where 



ttfe 



3 3^l-6Awiq' + p') ' 



(2 + ^)(l->g-V4-2g-2g^) / . ,A ^2^ 2^ 

«c = —, , , S = \/l - 6Aw[q +p> ■ 

6 (8 - 2^2 _ 4^4 _ 25 - 252 - 5V4 - 2S - 2S^) ^ 

(43) 

We now take the near-flat approximation by expanding to leading order in L = \Aw\ {(f' + 
p2) ^ 1. The mass, the horizon locations and the temperatures of the near-cold near-flat black 



hole are given by 

^cold 



rriQ 



I n 



out 

'o 



cosm 
I n 



V^+?(^l-j^+0{L')., 



.COM eextr , 525ao(g^ + p^)^ - 924ai(g^ + p^)p^ - 656a2P^ r <nfr^\ 



840(g2+p2)3 
105ao(g2 +p2)2(^_3 ^4^) ^ 1680ai(g2 +p2y^ ^ 1748a2P=^A 

840(g2+p2 



P - Pl=o + ^JKTZTT^^ L + U[L +A 



out _ out 105ao(g^+p^)^(3 + 4A) + 1680ai(g^+pVA + 1748a2P^A ^ , ^,^, , ^,, 
P - Pl=o 840(^2 + ^2)3 L + 0(L+A), 

s.extr J 1/2 / a2\ 

^on* ^ ™.t _ ll^A A , 5985ao(g^ +p^)=^ - 6090ai(g^ +p^y - 3464a2p^ \ 

^=° 87r v/^^T^ V ' 4620(g2 + p2)3 ) + 

+ 0(L2 + A2 + e2) , 

T^°^"^ = ^^'^' - 2L - L A2 / _ ^,^e.J. r 1/2 A + ^^ V 0(2,3/2 + ^4 ^ ^2) (44) 

47rv/g2+p2 \^ \ 2 J J 

where L = designates the near-extremal flat case results of the previous section. For A = 0, 
we get the cold near- flat limit. 



4 A Detailed Balanced Hofava- Maxwell Theory 

Detailed balanced Hofava gravity was presented in p^ . A detailed balanced version of the theory 
coupled to a scalar field, was derived in [4]. Here we discuss a detailed balanced version of 
Hofava gravity coupled to Maxwell and higher derivative magnetic terms. Detailed balance is a 
symmetry condition which serves to reduce the number of terms in the Lagrangian. This allows 
in our case to yield analytic non-perturbative solutions. 

The detailed balance Lagrangian density will be given by: C = Ck — Cy- The kinetic part is 
the same as before: 

j^^K = ^(^^.^^' - A^') + ::^9''{Ft^ - N'F,,){F,, - N'Fi,) . (45) 

10 



The potential part is defined to be of the form: 

1 . K^ ^ „..„., 1 



with the generahzed de Witt metric: 

Qijki = -{gikgji + gugjk) + _ ^. gtjgki , (47) 

and where 

for some three-dimensional Euclidean Lagrangian density W. For z = 3, W should include up 
to three spatial derivatives. 

The gravitational part of W, due to [1], is given by 

—=WHofava = 0^=^'^ 

Vg w^ ^fg 



X c- ' / 2 \ 

—pL^Uai-ava = - — —^""u ( djT f,^ + -T jj^'^km j + K^ " 2Avi/) , (49) 



where F* -^ are the Christoffel symbols. This is the Lagrangian density of topological massive 
gravity [23], with the first term being the gravitational Chern-Simons terno. Varying with 
respect to gij yields 

Using this, one gets the potential part in Cuofava dS]). 
For the electromagnetic part of W , we take 

-^Wem = ^^Ad.Ak + hF,,F'^ + 62^(AA, - R\Ai)d,Ak , (51) 

V 5* V 5' yg 

where 6i, 62 are constant parameters. This is the three-dimensional Chern-Simons-Maxwell- 
extended-Chern-Simons Lagrangian density in curved space [21]. Varying with respect to g^j 
yields 



^- ^ 7f ^ ^ "' ^'"'"'" - '*'^"^") " 



Akl 



+ b2—iF^"'VkFim + FimVkF'"' - gknFi^W'F^^) , (52) 



^Note that to be consistent with (jSO)) and (|H), we find this term to have an opposite sign compared to [T]. 
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where in the spherically symmetric case, the &2 term vanishes. Varying with respect to the Ai 
yields 

where in the spherically symmetric case, the hi and 62 terms vanish. Note that both variations 
yield gauge invariant expressions. Other terms may be added to Wem, which would yield gauge 
dependant expressions. Nevertheless, some of these may still lead to gauge invariant equations 
of motion. 

In the static spherically symmetric case, we substitute the ansatz f lT3|) . f lT^ into the La- 
grangian density to get 

^r^Creduced oc 3Awr' + 2fr' + 2fr^ - 2r^ + -^{2ffr'' - 2f'r'> - fr' + 2// - r^) + 
_ ^^'2^8 ^ 2p\^ + /^(2Awr^ - 2f'r' + 6fr^ - 6r^ - 5/3/) + 
+ ^(^K^rfr' + N^f'r' + N^fr') , (54) 

where /3 = fj^^^bi, we have set A = 1, and have omitted an overall factor independent or r. 

We now find non-perturbative solutions with higher order magnetic terms. For the projectable 
case with Nr ^ 0, N = 1, and constant /, we find the solution: 

AT ^^ l^w , , , ^ , 2m , (/ - 1)2 q2+p2 f 1,1-/, /3p' \ .,,. 

''^ = ^l^^'^^-'^-^^A^-^^^^' [-7^^A^^^2A^^) ■ ^''^ 

For the non-projectable case with Nr = 0, we find 

2 
A^' = / = -Awr^ + l + f3^- ^/-AAwmr + 2Awiq^+P^) ■ (56) 



These generalize the solutions with /3 = of pTl [T6] . 
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